Exact propagating nonlinear singular disturbances in strongly coupled dusty plasmas 
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The nonlinear longitudinal response ol a strongly coupled dusty plasma system is analytically 
investigated using the Generalized Hydrodynamic (GHD) model. It is shown that the Galilean 
invariant form of the model does not have soliton solutions, but can support a variety of nonlinear 
singular (yet conservative) solutions like shock structures of zero strength and propagating solutions 
with cusp like singularities. These novel entities should be detectable in numerical simulations and 
experiments studying large amplitude longitudinal excitations in such systems. 
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When micron sized dust is sprinkled in a plasma, the dust particles acquire a negative charge and shielding clouds 
large enough to balance thermal electron currents to ion thermal currents. The shielded dust particles surrounded by 
their ion and electron Debye clouds thus form an additional charged component of the plasma. If the dust density 
is large enough, the Yukawa interaction energy between two shielded dust particles may readily exceed their kinetic 
energies putting the dust component into a strongly coupled state whereas the electrons and ions still remain in a 
weakly coupled regime [H 0] . Strongly coupled dusty plasmas may be found in many physical situations ranging from 
planetary astrophysical plasmas Q to plasmas in industrial environment Q and even plasmas at the boundaries of 
thermonuclear fusion devices § . Such plasmas have also been produced in the laboratory for basic experiments where 
one has observed interesting phenomena such as formation and melting of crystalline structures [6|-|8| and excitation 
of collective modes in the crystalline as well as liquid phases of the dust component ■ There have also been a 

number of theoretical investigations of strongly coupled dusty plasmas using a variety of techniques [34I3- Among 
them the Generalized Hydrodynamic (GHD) model approach [16J, that augments the usual hydrodynamic description 
by a phenomenological visco-elastic term to account for the strong coupling correlations, has been fairly successful 
in describing the linear collective properties of a strongly coupled dusty plasma system. For example, the existence 
of transverse shear waves in the liquid state and strong coupling induced modifications in the dispersion properties 
of dust acoustic waves were predicted by the GHD model [16j and subsequently verified in laboratory experiments 
fl2l [iH ] and molecular dynamics simulation [l8[ . The application of the GHD model in the nonlinear high amplitude 
regime has however remained quite sparse and limited, primarily due to technical difficulties in analyzing the equa- 
tions. Unlike the weak coupling regime, the equations cannot be reduced to a KdV like paradigmatic equation that 
has soliton solutions. In fact, as we will demonstrate, the nonlinear GHD model cannot support soliton solutions. Our 
present investigations show, that under some realistic approximations and in the weak amplitude limit, the nonlinear 
GHD equations can be reduced to the Huntcr-Saxton equation, an intcgrable nonlinear PDE with a very different 
and distinct class of exact nonlinear solutions. In contrast to the usual KdV solitons, these solutions have a spatial 
singularity and yet retain conservative properties including many soliton like collision characteristics. In the more 
general case, for arbitrary amplitude excitations, we show that the GHD equations have nonlinear solutions in the 
form of propagating cusp like structures for the velocity perturbations. Our results therefore offer a new paradigm 
for analysing nonlinear excitations in strongly coupled dusty plasmas and point to new directions of research for 
numerical simulation and experimental explorations. 



The basic e quat ions of the Generalized Hydrodynamic (GHD) model describing the dusty plasma medium in one 
dimension are [16j . 
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Here, n is the dust density, u is the dust velocity along x and 77* is the viscosity coefficient. The evolution of u 
incorporates characteristic traits associated with strongly coupled matter through a relaxation parameter r m . For 
time scales longer than r m the medium behaves like a viscous liquid, whereas at shorter time scales the memory effect 
persists and the system shows solid like elastic properties. The operator d/dt = d/dt + ud/dx denotes the total time 
derivative. It may be mentioned here that the convective term appearing with coefficient r TO in Eq.([2]) has often been 
dropped as a simplification in several earlier studies [Ul23|- However, this term ensures the Gallilean invariance of 
the equation and ought to be retained. The scalar potential <p is determined from the Poisson's equation (Eq.Q). 
The time scales associated with the evolution of the dust being long, the electron and ion species are assumed to 
satisfy the Boltzman distribution n e — [i e exp (<7j0) and rij = p-iexp (—</>) and have been accordingly specified in the 
Poisson's equation. The dust pressure is modelled by a simple equation of state P = ksnTd. Above equations have 
been written in terms of normalized fields described in one of our earlier works [2l| . 

In this paper we concentrate on the strongly coupled limit k^/r m rf » 1, where k is the inverse scale length of the 
solution. We shall primarily use analytic methods for drawing conclusions. 

If ff /r m >> C\ a (Cda being the dust acoustic speed), i.e. when the elastic wave dominates the dust acoustic speed, 
one can ignore the contribution from the scalar potential </> and the thermal contribution due to P in the momentum 
equation. Physically, this is the regime when clastic coefficients due to correlations dominate over Boltzmann screening 
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and thermal dispersion effects. The dusty plasma medium is, however, still in a fluid molten state with no lattice 
formation. In this limit the dust fluid is governed by the following simplified equation: 
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It should be emphasized that compressional velocity perturbations in the dust fluid will still produce density distur- 
bances, which in turn will be shielded by electrons and ions producing potential perturbations. The inequality at the 
beginning of this paragraph ensures that the reaction back of these driven disturbances on the momentum equation is 
negligible. Physically eq.Q contains dispersion free linear elastic waves that are supported by the correlation driven 
elasticity coefficient and nonlinear contributions through inertial effects appear through the convective terms. In 



principle, linear wave dispersion may be introduced through a k dependent form of r m 22(; here we assume that this 
effect is small. Note that the second convective derivative which arose through constraints of Galilean invariance is 
playing a crucial role in the nonlinear dynamics. This equation can also model plastic flow deformation disturbances 
in solids undergoing failure through severe stresses. 

In the weakly nonlinear regime Eq.ffl can be subjected to a reductive perturbation analysis, by expanding, u = 
A + euS 1 ' + e 2 u^ + . . . and v = e 3 ^ 2 (v^ + ev^ + . . .) where v = (d/dt + ud/dx)u. Further, using the stretched 
variables, £ = e 1 / 2 [x — \t), r = e 3 / 2 i, taking A = \Jrf jr m , and retaining terms upto second order for the u and v 
fields, we can obtain the following single equation in the variable u^ (rewritten below without the superscript), 

(u T + uu S ) 6 = -u\ (5) 

The left hand side equated to zero is the nonlinear equation for dispcrsionless waves with the convective nonlinearity 
giving indefinite steepening of waves which can lead to wave breaking or form shocks, or solitons depending on whether 
nonlinearity, viscous dissipation (Burgers equation) or dispersion (Korteweg de Vries equation) dominates the physics 
of steepened waves. Here the extra convective derivative nonlinearity of the simplified generalized hydrodynamic 
model equation (£[]) is responsible for the nonlinear term on the right side of equation (|5|). This term dramaticaly 
changes the character of the equation and the nature of its solutions. 

Equation ([5]) is the so-called Hunter-Saxton equation, which has been derived earlier [23] for director fields in 
liquid crystals, where the positional disorder of polymer molecules gives the medium fluid properties whereas the 
orientational order due to correlations gives them crystal like properties. It is also the high frequency limit of the 
Cammasa- Holm equation (2~H , which has been derived to describe the nonlinear dynamics and wave breaking of 
shallow water waves. These equations belong to a new class of equations which can be derived from variational 
principles in more than one non equivalent forms. They typically have an infinite number of conservation laws and 
possess singular solutions with infinite derivatives. If these solutions are propagating, they pass through each other 
undisturbed, except for a phase shift, somewhat like solitons. 

We now recapitulate some properties [25j j of the Hunter Saxton equation and its solution, which are of relevance to 
our problem, Firstly, integrating Eqn. §5§ over £ we note that because of the positive definite value of the integral 
on the right side, if the solutions leave one boundary unperturbed, the other boundary is perturbed, thus showing 
the impossibility of smooth periodic or isolated solutions with undisturbed boundaries. Secondly, a Lagrange variable 
treatment for the derivative u^ shows that the velocity derivative blows up in a finite time. This indicates that the 
nonlinear disturbances will lead to a wave breaking like behaviour in a finite time. However, since correlations lead 
to elasticity, we find wave breaking phenomenon with a difference. This is best illustrated by the exact solution (23j 
below. 

Hunter Saxton (HS) equation has step like piecewise continuous and non - smooth solutions. A single step solution 
may be written as 

u(£,T;a,P) =U(£,r,a) if r<0 

= U(£,t;P) if r>0 (6) 

where a and f3 are positive constants with condition /3 < a and, 
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U(£, t; a) = —ar — oo < £ < 

= < £ < oo (7) 



4 



l I V Tj <T 2 <0<T 3 <T 4 




P = tx icunservuiivi' 





jS = o, cc ^ | dissipati ve solution 



FIG. 1. Schematic view of time evolution of solution ([6]) of the Hunter Saxton equation for the conservative and dissipative 
cases (redrawn from HU). 

Also £/(£ , 0, ; a) — so that U is a continuous function of £ and r. Mathematically, the solution is a weak solution of 
the HS equation satisfying the condition 



is dissipative but still weakly admissible. The solutions are non propagating in the wave frame travelling to the right 
with the linear phase velocity y/ (r?*/r m ). 

Unlike the inviscid Burger equation, where the step becomes vertical (shock solution) and acquires a finite steady 
value consistent with conservation laws, the HS equation has a step solution with a slope and a step size that are 
time dependent. The remarkable feature of the solution is that as the left and right corners of the linear segment 
having negative slope collide for the creation of a 'shock wave' with infinite slope (which may lead to wave breaking) 
the spatial support in real space diminishes to a point and the step size vanishes simultaneously. Since the region 
of transition diminishes to zero size as the step approaches verticality, J u^dx can remain conserved. In this case no 
norm is lost and 'energy' is conserved. This is unlike normal Burger's like shock wave where the step is constant and 
some energy is converted to heat. In fact at complete verticality, the HS solution has no step and that is why it is 
sometimes called a 'shock wave of zero strength'. However, we note that such shock waves of zero strength can form 
cither conservative (/3 = a) or dissipative (/3 < a including j3 = 0) global solutions. For example, for a ^ 0, /3 = 
the disturbance starts with a positive step of u on the left at r < ( see Fig. [I} and then this step goes to zero 
at t — when the step becomes vertical. Thereafter (r > 0), the disturbance vanishes from everywhere. This is a 
dissipative global solution for which the entire energy in the initial disturbance damps away and disappears. This is 
akin to conventional wave breaking, where the infinite slope leads to toppling of the wave and conversion of coherent 
wave energy into chaotic multi stream motions. The conservative global solution, on the other hand, corresponds to 
j3 = a and results in a fresh disturbance with positive slope at r > 0, where the conserved J u^dx energy results in a 
diminishing slope disturbance in a widening region (Fig.[TJ) as r increases . This is a remarkable sequel to wave breaking 
with infinite slope at t = 0, a sequel in which the entire elastic energy of the infinite slope wave trapped in a region 
of zero size reappears as a coherent elastic disturbance. This is only possible due to the presence of the RHS of HS 
equation. The HS equation thus shows evidence for self consistent nonlinear elastic waves supported by correlations, 
which steepen indefinitely but peter out in strength before they reach verticality; these are weakly nonlinear waves 




(8) 



for arbitrary test function $; this is weakly admissible if 
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which steepen and want to break but cannot do so because of strong coupling and correlations. Physically, one may 
picture the global dissipative and conservative solutions described above as longitudinal disturbances along a spring 
attached to a wall. Imagine a compressional disturbance coming towards the wall, steepening and becoming infinitely 
compressed at the wall. The subsequent behaviour can be either inelastic (global dissipative solution) with the entire 
energy in the disturbance dissipated at the wall (say, because of plastic failure of the spring) and nothing returning 
back or elastic with a longitudinal disturbance of equal magnitude returning from the wall. 

We now examine the arbitrary amplitude Eqn.Q, for some special exact solutions. We wish to explore the possibility 
of singular solutions with singular derivatives. For this purpose we introduce a functional of x — j3t as the independent 
variable, viz. 

f = k(x - f3t) + G(f) (10) 

and assume that u is a function of x, t only through the combination /. This means solution u(x, t) is time independent 
in a frame moving with velocity (3. If G(f) is chosen as a constant, it simply acts as a phase term. If G, on the other 
hand is chosen as a function of /, then / becomes a functional of x — j3t; this dependence can be made implicit if G 
depends on / through u. This is indeed what we shall do, so as to derive an exact cusp like solution of u{x, t) with 
infinite derivatives at the maximum. From Eqn. (|10p we get, 



d u — f3 du rj* d 1 du 
{U ~ P) ~df 1 - Gj ~df = ~^df l-Gf'aJ (U) 



Introducing u — /3 = y/rj*/T m U and making the choice 1 — Gf = ^/rf /T m U, we get the equation 

TT d 2 U _ d 1 dU 

u W~dfuW (12) 

This equation can be integrated once and rewritten as an energy integral 

where the constant of integration has been absorbed in function / as a multiplication factor and sets the scale of the 
solution. Equation (| 13[) is the equation of a zero total energy effective particle moving in a potential energy bowl 
with an inverse parabolic form near the origin that blows up at U = ±1. Thus if the effective particle starts from 
U = with near zero velocity (W), it falls towards U = 1 slowly at first and then rapidly, till it reaches the "wall" 
(singularity at U = 1) from where it is reflected back . It takes an infinite time to climb back to U = again. The 
resulting solution is an isolated soliton like solution with a cusp at the maximum U = 1. The solution showing this 
property is 



U = sech(f + y/l - U 2 ) (14) 

This solution is illustrated in Fig [5] and shows infinite derivatives and a cusp singularity at U = 1 (that is, 
u = f3 + (rj* /r m ). The phase speed and the scale size of theese solutions are not directly determined by the 
maximum amplitude. The nature of singularity can be explored by expanding the solution around U = 1 and 
demonstrating that U' ~ Z" 1 ^ 3 - These cuspon like solutions are steady propagating solutions which are singular 
at a point and are dithering at the wave breaking amplitude. Physically, such solutions might arise when smooth 
nonlinear waves acquire amplitudes close to wave breaking, but because of conservation laws squeeze the infinite 
derivative region to a point with a finite elastic energy content. That this is indeed so, can be ascertained from the 
integral / U' 2 df = J U'dU = J UdU j\J\ — U 2 = 1 in normalized units. Such finite energy content singular solutions 
may have special stability properties. It may be worthwhile to look for an infinite number of conservation laws for the 
exact Eqn.Q, for if they are found then two such solutions could propagate through each other without distortions, 
like solitons. In fact, the stability and accessibility of these solutions from arbitrary initial conditions is a topic worthy 
of further investigation. 

In conclusion, we have investigated the nonlinear longitudinal disturbances in a strongly coupled dusty plasma 
using the description of generalized hydrodynamic equations. We work in the limit of very strong coupling where the 
correlation induced elasticity dominates the dispersion relation for linear longitudinal waves. For the weakly nonlinear 
case we have found that the nonlinear waves are described by the Hunter Saxton equation. This equation has exact 
non smooth solutions which are in the nature of time dependent steps and slopes propagating at the speed of linear 
clastic waves. We have considered a single step solution and shown how it gives rise to a shock wave of zero strength, 
which either produces a global dissipative solution or a global conservative solution. The former leads to a damping 
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FIG. 2. Form of potential energy "bowl" and zero energy "particle" orbit corresponding to a cuspon solution. 



of the coherent wave energy like in normal wave breaking events whereas the latter leads to the reappearance of a 
step like elastic disturbance which carries off the energy trapped in the zero size singularity region associated with 
the shock wave of zero strength. For the arbitrary amplitude case we have discovered the existence of cuspon like 
solutions, viz. propagating longitudinal isolated disturbances with cusp like singularities at the maximum. It would 
be very worth while to look for cuspons, shock waves of zero strength and global dissipative and conservative solutions 
in laboratory and simulation investigations of nonlinear disturbances in strongly coupled dusty plasmas. It is possible 
that the cuspon solutions are related to cusp like solutions already observed in some dusty plasma experiments [26j . 
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